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The use of fused-silica ribbons as suspensions in gravitational wave inter- 
ferometers can result in significant improvements in pendulum mode thermal 
noise. Surface loss sets a lower bound to the level of noise achievable, at what 
level depends on the dissipation depth and other physical parameters. For 
LIGO II, the high breaking strength of pristine fused silica filaments, the cor- 
rect choice of ribbon aspect ratio (to minimize thermoelastic damping), and 
low dissipation depth combined with the other achievable parameters can re- 
duce the pendulum mode thermal noise in a ribbon suspension well below the 
radiation pressure noise. Despite producing higher levels of pendulum mode 
thermal noise, cylindrical fiber suspensions provide an acceptable alternative 
for LIGO II, should unforeseen problems with ribbon suspensions arise. 
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I. INTRODUCTION 



One of the most important limitations to the sensitivity of long baseline gravitational 
wave detectors such as LIGO [|TJ, VIRGO 0, GEO 600 §,and TAMA [g is thermal noise 
associated with the test masses and their suspensions. Designs for advanced detectors pro- 
pose either fused silica or sapphire test-masses. For fused silica test-masses internal mode 
thermal noise is expected to be an important source of noise from approximately 20 Hz to 
a few hundred Hertz, whereas pendulum mode thermal noise is more important below this 
range. 

Pendulum mode thermal noise is due primarily to dissipation in the suspending filaments. 
It is imperative, therefore, to minimize the intrinsic losses in the filament. Many current 
detector designs use metal wires to suspend the test-mass, but metals, with > 10~ 5 , result 
in unacceptably high levels of pendulum mode thermal noise. Fused silica has lower loss 
((f) pa 3 x 10~ 8 ) and monolithic fused silica suspensions have been shown to have much 
higher Q than metal wire suspensions Such a monolithic suspension system is being 

developed and adopted for use in the GEO 600 detector 0, while variations of this design 
are being considered for LIGO II. In particular, fibers with circular cross sections may be 
replaced with fused silica ribbons [|3|-|T0|j, which allow the suspension filaments to be very 
thin and compliant in the direction of motion (|TT|-[T3|. 

Experiments indicate that in thin fused silica filaments much of the dissipation takes place 
in a layer near the surface [0. The level at which surface loss affects the total pendulum 
mode dissipation depends on the filament thickness and geometry, and influences the choice 
of suspension design parameters. To investigate this, we calculate the pendulum mode 
thermal noise, including surface dependent loss, as a function of the design parameters for 
fibers and for ribbons. 



In the absence of external sources of dissipation, the dissipation in the fundamental mode 
of a pendulum suspended by a filament is given by JT^| 



where is the loss angle of the unloaded suspension filament, Y is the Young's modulus of 
the filament material, I is the cross-sectional moment of inertia, M is the supported mass, 
g is the acceleration due to gravity, and L is the filament length. In general, <ft an d hence 
$ will be a function of frequency. The coefficient of is called the dissipation dilution 
factor and is the ratio of elastic energy (subject to dissipation) to the total energy stored 
in the pendulum mode, which is predominantly gravitational energy (non-dissipative). The 
right hand side of this equation should be multiplied by two if the mass is constrained from 
rotating in the plane of oscillation, since bending then occurs both in the region where the 
filament leaves the support and in the region where the filament leaves the mass. If the test 
mass is suspended by N filaments Eq. |l] should be multiplied by y/N. For fibers of diameter 
df, and ribbons of thickness d r and width w, we have 



II. DISSIPATION DILUTION 
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where subscripts / refer to fibers and subscripts r to ribbons. Rewriting Eq. |l] using these ex- 
pressions for the cross-sectional moment of inertia, allowing multiple filaments and assuming 
that the suspension constrains the filaments to bend at both ends, we have 



$=<( v ; 4M ^ 0/ ' (3) 

Tm0 & ribbons, 

where M is the suspended mass, Nf is the number of suspension fibers, 0/ is the loss angle 
of the unloaded fibers, N r is the number of suspension ribbons, and <p r is the loss angle of 
the unloaded ribbons. The limit to how much dissipation dilution we can obtain, and hence 
the lower limit to the pendulum mode thermal noise, is set by the values obtainable for the 
parameters in these equations. They are limited by a number of material and technological 
concerns, especially by the value achievable for the loss angle <pf or <p r . This loss angle may 
depend on a number of factors including the bulk material loss angle, surface loss, and the 
filament geometry. 

If <pf and <p r were independent of filament thickness and roughly equal, Eq. [3] would 
indicate that by using very thin but wide ribbons one could obtain lower dissipation $, and 
hence less pendulum mode thermal noise, than by using fibers of similar load bearing ca- 
pacity. However, since surface loss becomes increasingly important as the filament thickness 
is reduced, the enhanced dissipation dilution obtainable using thin ribbons is moderated by 
an increase in <b r . 



III. THERMAL NOISE IN THE PRESENCE OF SURFACE LOSS 

The loss angle for a sample, including surface loss, may be expressed as 

= 0bulk(l + / i 7/^)> ( 4 ) 

where 0tmik is the loss angle of the bulk material, // is a geometrical factor and d s is the 
dissipation depth which parametrizes the filament size at which surface loss becomes im- 
portant. The geometrical factor /i describes the emphasis placed on the condition of the 
surface due to the sample geometry and mode of oscillation while the dissipation depth d s 
describes the amount of surface damage and the depth to which it penetrates. Equation |] 
serves to define d s , whose value for a given sample may be determined by experiment. The 
geometrical factor is given by 

_V j s e\r)d\ 
* S / v e 2 (f)dV K ' 

where r denotes a point in the sample, e(r) the strain amplitude, V is the volume of the 
sample, S the surface area of the sample, V is the set of points comprising the volume, and 
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S is the set of points comprising the outer surface. For transverse oscillations of fibers and 
ribbons we have 



(3 + a)/(l + a) 



fibers, 
ribbons 



(6) 



where a is the aspect ratio of the combined ribbons, a = d r /W with W = N r w being the 
total combined width of the ribbons. 

Experiments suggest that bu ik is approximately constant over the frequency range of 
interest for LIGO @,^,^]. For simplicity, we will assume 0buik to be constant. Substituting 
Eqs. f| and || into Eq. |3] we obtain 



Y/16aL 2 (df + 8<i s )0buik fibers, 
Y fYloL 2 id r + (6 + 2a)<i s )0buik ribbons, 



(7) 



where a is the filament stress. In both cases, the first term is the traditional expression for 
dissipation dilution, while the term involving d s represents a reduction of the dilution due 
to the increasing importance of surface loss as the filament thickness is decreased. For very 
thin filaments the term involving d s dominates and the loss angle becomes independent of 
the filament thickness. 

From the fluctuation-dissipation theorem ||15|| , we find the power spectrum of the pendu- 
lum mode thermal fluctuations above the pendulum mode resonance, at angular frequency 
to: 



2 4:k B Tg 



(8) 



where 10 <; J g/L' } T is the temperature of the suspending filaments and V is the radius of the 
arc traced out by the center of mass during pendulum mode oscillation. For convenience, we 
will take V « L. Inserting $ from Eq. [7|, and including the contribution from thermoelastic 
damping we have the expression for the pendulum mode thermal noise: 
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The thermoelastic damping term t h is given by |16 
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(10) 



where a is the thermal expansion coefficient of the filament material, C is the heat capacity 
per unit volume, and D is the thermal diffusion coefficient. 



3 



IV. ADVANCED INTERFEROMETERS 



Using Eq. |S], we can now make estimates for the level of pendulum mode thermal noise 
achievable in advanced interferometers and investigate the dependence on filament thickness 
and geometry. Using the results of previous experiments and design studies, most of the 
parameters can be well bounded with some reasonable assumptions. From these parameters, 
we can obtain upper and lower bounds on the pendulum mode thermal noise at a given 
frequency as a function of ribbon thickness. This analysis assumes that losses extrinsic to 
the filaments (e.g. recoil of the suspending structure or lossy filament-to-test-mass bonds) 
have been made negligible. 

The achievable fiber diameter depends on the achievable stress a and on the mass M. 
The fiber diameter in Eq. R| can be replaced by 



d f = ^4Mg/nN f a. (11) 

The remaining parameters M,L,Nf, a, and d s are independent and the achievable pendulum 
mode thermal noise depends on the bounds established for these parameters. 

It is clear from Eq. ||] that an efficient way of minimizing the thermal noise is to make the 
length of the suspension as large as possible. However, the value of L is bounded above by 
the minimum allowable spacing f min of the violin mode frequencies. The frequency spacing 
must be kept above about 300 Hz to allow reasonably large intervals of the spectrum to be 
free of violin modes. The frequency spacing limited L is 

(12) 

where / m j n is the minimum allowable spacing of the violin modes. The range of possible 
lengths is determined by the range of possible stress to which the filaments will be subject. 
This in turn depends on the breaking strengths achievable for fused silica filaments. Many 
measurements have been reported on the breaking strength of fibers manufactured from nat- 
urally occurring, and synthetic, vitreous silica |TT|,|T8f. Little is known about the strength of 




ribbons, though one is tempted to assume their strengths are similar. Values reported for 
the breaking strength of fibers in tension at room temperature vary greatly depending on the 



condition of the fibers ||19|| , but strengths on the order of several gigapascals at room tem- 
perature, in fibers with diameters as large as 1 mm have been reported |2"0"|--f2"!^| . By assuming 
that the filaments are only loaded to a fraction of their breaking strength we assign the range 
of possible stress to which the filaments will be subject as 0.1 GPa < a < 1.0 GPa. Substi- 
tuting these values into Eq. [12] we obtain the range of possible lengths 0.36 m < L < 1.1 m. 
In principle, the physical design of the suspension also places an upper limit on the length, 
but ultimately this limit is likely to be less stringent than that due to the frequency spacing. 

For the number of filaments we will choose Nf = N r = 4. This reflects the most likely 
choice for the suspensions in advanced detectors which require "marionette" control of the 



test masses |10|. Analytically, the number of ribbons does not enter into the calculation as, 



for a given stress a and ribbon thickness d r , only the total combined width of the ribbons 
W is fixed. 
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In order to avoid excessive radiation pressure noise in LIGO II, the suspended test masses 
must have a mass M of 30 kg. However, if this is not feasible, the LIGO I mass of 10 kg can 
be used as a fall-back. We take 10 kg < M < 100 kg to allow for possible advanced designs 
that utilize even larger masses |10| . 



For the bulk material dissipation in fused silica, Gretarsson and Harry have measured 
0buik = 3.3 ± 0.3 x 10~ 8 [[/J. Others have measured similar values for the bulk material 
dissipation in samples of different geometry P,[Z5|-|2"5| . We shall adopt the relatively reliable 
upper limit 0buik = 3.6 x 10~ 8 and somewhat uncertainly set the lower limit as </>buik = 
2.5 x 10~ 8 

Finally, for untreated fused silica fibers drawn in a natural gas flame, d s = 180 ± 20 /im 
has been measured fTj. It should be noted that the factors resulting in a given quality 
of fiber surface are not well quantified. Fibers pulled from silica rods with different initial 
surface conditions, or fibers drawn using a different production method, may have a different 
dissipation depth or surface loss than that found in the measurements above. The geometry 
of a filament could also have some effect on the quality of the surface layer, e.g. through 
different cooling stresses during fabrication, and our assumption that fused silica ribbons 
have the same surface properties as fused silica fibers has not been tested. However, given 
these assumptions we set an upper bound of d s = 200 /mi, which should be reliable. To 
estimate a lower bound for d s we will use a Q-measurement of a ribbon of thickness 50 /im, 
made of natural fused quartz ||. One mode of this ribbon showed a Q much higher than 
others. After subtracting the loss due to thermoelastic damping (Eq. [H]) and assuming the 
remaining loss is mainly surface loss, the equivalent d s for a similarly limited fused silica 
ribbon can be estimated at 30 /im. We, therefore, set the range of possible dissipation depths 
at 30 /im < d s < 200 /im. 

Table | summarizes the best and worst case estimates for the parameters, and also a "best 
guess" for the most probable values. Figure [I] shows the levels of pendulum mode thermal 
noise at 10 Hz as a function of filament thickness for each of the three sets of estimates of the 
parameters. The graphs for ribbon filaments all show a maximum around 400 /im. This is 
the thermoelastic damping peak, which for all but the most optimistic case must be avoided 
if the desired levels of pendulum mode thermal noise are to be achieved. Figure [1] also shows 
the value of using ribbons rather than fibers as suspension filaments. Since the diameter 
of fibers cannot be independently reduced, the pendulum mode thermal noise is dominated 
by thermoelastic damping. The use of ribbons allows us to evade this problem by residing 
in the surface loss limited regime. To evade the thermoelastic regime and to obtain better 
dissipation dilution one might be tempted to use the thinnest ribbons possible. However, at 
small thicknesses, below the thermoelastic peak, the graphs begin to level off. This is due 
to surface loss which sets a minimum to achievable pendulum mode thermal noise of 

2 2Ak B Tg [y~ 

It is clear from the plots that, for all but the most optimistic values of d s , reducing the ribbon 
thickness below about 50 /tm (corresponding to individual ribbon widths w of 5 mm, 3 mm, 
and 5 mm for the three cases) does not result in significant reductions of pendulum mode 
thermal noise. To satisfy LIGO II requirements, the pendulum mode thermal noise at 10 Hz 
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should be less than about 2 x 10 m/yHz [ 1Q| . Even in the presence of surface loss, only 



the worst case scenario does not achieve this level. In the most probable case, pendulum 
mode thermal noise will be lower than other noise sources at 10 Hz (radiation pressure 
noise, fused silica internal mode thermal noise), provided the ribbon thickness is kept below 
the thermoelastic regime. In the most optimistic case, other noise sources dominate the 
total noise at 10 Hz regardless of ribbon thickness. The most probable estimate for fiber 
suspensions gives pendulum mode thermal noise that is just acceptable for LIGO II. If 
there are unforeseen problems with ribbon suspensions, fiber suspensions may still prove an 
acceptable alternative. We reiterate that in the comparison of fibers and ribbons, we have 
assumed that the breaking strength of fibers is not significantly greater than that of ribbons 
of equal cross-sectional area; we have also assigned them identical surface properties. Further 
research is required to test these assumptions. Research is continuing on ribbon suspensions 
within the LIGO research community, and additional emphasis on surface properties and 
breaking strength is warranted. 



V. COMPARISON OF LOW FREQUENCY NOISE SOURCES IN LIGO II 

While studying the pendulum mode thermal noise at 10 Hz is a good way to gain insight 
into the effect of the different physical parameters on the level of this source of noise, 
the comparison with other sources of noise must be done over the entire range of relevant 
frequencies. 

For clarity, we now specialize to a single set of values for the physical parameters of 
the suspension filaments-those proposed for LIGO II [ rP[ . These parameters are shown in 



the last column of Table |. With a four ribbon suspension, it follows from the proposed 
stress that each ribbon should have a cross-sectional area of 5.5 x 10~ 7 m 2 , giving a width 
of 5.5 mm for the 100 /im ribbon thickness proposed. In general, the values for all these 
parameters fall between the worst case and most probable case scenarios. As such, the LIGO 
II proposal is fairly conservative, and better noise performance may be achieved. 

The LIGO II proposal does not, however, specify a value for d s . In keeping with the 
conservative spirit of the other parameters, we choose d s = 180 /xm. This falls between the 
worst case and most probable case and corresponds to the measured value in fibers without 
strict handling requirements 0. 

Figure shows the pendulum mode thermal noise of a fiber suspension in relation to 
estimates for the other sources of noise in the interferometer. Figure |3] shows the same 



comparison for the pendulum mode thermal noise of a ribbon suspension |[26|| . Note that 
the noise due to radiation pressure in the LIGO II design is greater than the pendulum 
mode thermal noise for the ribbon suspension. For greater low frequency sensitivity, radi- 
ation pressure can always be reduced by lowering the amount of laser power at the beam 
splitter. This could reduce the noise to the pendulum mode thermal noise level in the low 
frequency band but will increase the amount of laser shot noise at higher frequencies. From 
figures |2| and |3] it is clear that while a ribbon suspension leads to lower pendulum mode 
thermal noise than a fiber suspension, the pendulum mode thermal noise for the fiber sus- 
pension is still comparable to radiation pressure noise in the relevant frequency band. If 
there are unforeseen problems with ribbons (buckling, lower strength, etc.), fibers do provide 



6 



an acceptable, if less attractive, alternative. 
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TABLES 



worst case 



most probable 



best case 



LIGO II proposal 



0bulk 

d s 
a 
M 
L 



3.6 x 10" 
200 

0.1 GPa 
10 kg 
0.35 m 



3.3 x 10- 
100 fjan 
0.5 GPa 
30 kg 
0.79 m 



2.5 x 10~ 8 
30 fim 

1.0 GPa 
100 kg 

1.1 m 



3.3 x 10-* 
180 jum* 
0.13 GPa 
30 kg 
0.39 m 



TABLE I. Estimates for the parameters determining the pendulum mode thermal noise, 
"not specified in proposal 
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FIG. 1. Solid lines: Worst case (column 1 in Table |), most probable (column 2 in Table |) and 
best case (column 3 in Table |) estimates for the level of pendulum mode thermal noise for ribbons 
vs. ribbon thickness. Dashed lines: Same, but neglecting surface loss. Asterisks: Worst case 
(column 1 in Table |), most probable (column 2 in Table |) and best case (column 3 in Table |) 
estimates for the level of pendulum mode thermal noise in 4-fiber suspensions with maximally 
loaded fibers. In the case of ribbons, the minimum thickness is set by requiring the total combined 
width of the ribbons W be less than or equal to 3 cm. The maximum ribbon thickness is set by 
requiring d r /W < 0.5. 
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FIG. 2. Frequency domain comparison of pendulum mode thermal noise of a fiber suspension 
with other low frequency noise sources in LIGO II. The assumptions in the LIGO II white paper 
concerning the output optics have been slightly modified to reflect the addition of thermoelastic 
damping in the test masses. The noise is expressed as motion of the face of a single test mass. 
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FIG. 3. Frequency domain comparison of pendulum mode thermal noise of a ribbon suspension 
with other low frequency noise sources in LIGO II. The assumptions in the LIGO II white paper 
concerning the output optics have been slightly modified to reflect the addition of thermoelastic 
damping in the test masses. The noise is expressed as motion of the face of a single test mass. 
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